Symmetry-allowed phase transitions realized by the two-dimensional fully frustrated 

XY class 
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A 2D Fully Frustrated XY(FFXY) class of models is shown to contain a new groundstate in 
addition to the checkerboard groundstates of the standard 2D FFXY model. The spin configuration 
of this additional groundstate is obtained. Associated with this groundstate there are additional 
phase transitions. An order parameter accounting for these new transitions is proposed. The 
transitions associated with the new order parameter are suggested to be similar to a 2D liquid- 
gas transition which implies ^-Ising like transitions. This suggests that the class of 2D FFXY 
models belongs within a U(\) <g> Zi <8> /^-designation of possible transitions, which implies that there 
are seven different possible single and combined transitions. MC-simulations for the generalized 
fully frustrated XY (GFFXY) model on a square lattice are used to investigate which of these 
possibilities can be realized in practice: five of the seven are encountered. Four critical points are 
deduced from the MC-simulations, three consistent with central charge c = 3/2 and one with c — 1. 
The implications for the standard 2D FFXY-model are discussed in particular with respect to the 
long standing controversy concerning the characteristics of its phase transitions. 



I. INTRODUCTION 

The two-dimensional(2D) fully frustrated XY 
(FFXY)-model describes a 2D Josephson junction 
array in a perpendicular magnetic field with the strength 
of the magnetic field corresponding to one magnetic flux 
quanta for every second plaquette of the array. The 
phase transitions of this model on a square lattice have 
been the subject of a long controversy^ - — .The emerging 
canonical picture is that the model has two relevant 
phase ordering symmetries: an angular £Y(l)-symmctry 
and a ^2-chirality symmetry^^*ii. As a consequence, 
the model has often been assumed to belong within 
the designation 14(1) ® Z^&^'-H. The controversial 
questions have been: Does the model undergo a single 
combined transition or two separate transitions and, if 
the latter, in which order do the transitions occur? The 
emerging consensus is two separate transitions: as the 
temperature is increased first a Kosterlitz-Thouless(KT) 
transition associated with the angular U(l)- symmetry 
and then at a slightly higher temperature a Z2-chirality 
transition^. The cause of the controversy can, ret- 
rospectively, be attributed to the fact that the two 
transitions are extremely close in temperature. 

We here generalize the 2D FFXY model into a a wider 
2D FFXY-class of models by changing the nearest neigh- 
bor interaction in such a way as to keep all symme- 
tries. This generalized 2D FFXY-class is shown to con- 
tain an additional groundstate. The existence of this ad- 
ditional grounstate leads to a phase diagram containing 
four sectors.—. We here show that it has seven different 
phase transitions lines and four multicritical points. We 
use Monte Carlo simulations to establish the characters 
of the transitions of this phase diagram. Our simulations 
suggest that three of the critical points are consistent 
with the central charge c = 3/2 and one with c = 1. 

In section 2 we define the 2D FFXY-model and in sec- 



tion 3 we describe the structure of the new ground state. 
In section 4 we propose an order parameter associated 
with the transition into this new groundstate. In section 
5 we give the results for the various phase transitions ob- 
tained from Monte Carlo simulations and determine the 
character of the four multicritical points by invoking a 
relation between the central charge c and the bulk criti- 
cal indices. In section 6 we discuss the original 2D FFXY 
model in view of our results. We also comment on related 
models not contained within the class of fully frustrated 
XY model discussed in the present investigation. Finally, 
some concluding remarks are given in section 7. 



II. GENERALIZED FULLY FRUSTRATED XY 
MODEL 

The Hamiltonian which defines the 2D fully frustrated 
XY-class models on an L x L square lattice is given by 



H = ^2u(^ ij =9i-e j -A 4j ), 



(1) 



with 4>ij S [ —7r i , where the sum is over nearest neigh- 
bor pairs. The phase angle 9i for the ith site at the lattice 
point (xi, yi) satisfies the periodicity 6i + L± — &i+Ly = 
The magnetic bond angle Ay is defined as the line inte- 
gral along the link from i to j, i.e. Ay = (27r/<I>o) f J A-dl 
with the magnetic vector potential A for the uniform 
magnetic field B = Bqz in the z direction. With the 
Landau gauge taken, Ay = 2nfxi for the vertical link 
and Aij = for the horizontal one, where the frustra- 
tion parameter / measures the average number of flux 
quanta per plaquette. The fully frustrated case corre- 
sponds to / — 1/2 with a half flux quantum per plaquette 
on the average. The Boltzmann factor, which determines 
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FIG. 1: Interaction potentials (7(0) in Eq. 2 at various values 
of p are compared in (a). The standard XY model correspond- 
ing to p = 1 is also compared with the Villain interaction po- 
tential in (b). All interactions have the same symmetry and 
have the identical quadratic form at small <j>. 



the thermodynamic properties, is given by exp(— H/T) 
where T is the temperature. The interaction potential 
U(<fi) — U(4> ± 2tt) is periodic in 2ir and is quadratic to 
lowest order in cf> so that U(4>) ~ <\> 2 . These conditions 
for the interaction potential defines the class: the mem- 
bers of this class are distinguished by the explicit form of 
the interaction potential U (4>) . If the relevant symmetry 
class is U(l) <g> Z2, then in principle three transitions are 
possible: separate IA(V) and Z2 transitions or a merged 
14(1) (g) Z 2 transition. However, the number of allowed 
phase transitions for the FFXY-class is much larger 13 . 
The implication is that by just changing the specific form 
of U((j)) within the FFXY-class one could encounter a 
plethora of phase transitions. In order to verify this, 
we choose a parametrization of U(<p) and find the phase 
transitions corresponding to this parametrization using 
Monte Carlo simulations techniques. This strategy was 
employed earlier in Rcf. 13 . The parametrization is of the 
form U(4>) wher e 14,15 

C^) = l[l-cos 2 ^)] (2) 
P ^ 

and p = 1 corresponds to the standard FFXY since 
2[1 - cos 2 (0/2))] = 1 — cos(</>). The members of the 
FFXY class, which belong to this parametrization, was 
in Refii 3 - termed the Generalized Fully Frustrated XY 
(GFFXY) model. Figure la shows a sequence of interac- 
tion potentials U((f>). 

To sum up: The 2D FFXY class which we discuss here 
is obtained from the standard 2D FFXY by generalizing 
the interaction potential within the allowed conditions: 
U(cj)) is a monotonously increasing function in the inter- 
val <f> S [0, 7r], U(4>) — U(4> ± 27r) is periodic in 27r and 
is quadratic to lowest order in <fi so that U(4>) ~ <fi 2 ■ 
The GFFXY model is by construction contained within 
this class. The Villain interaction is also contained in 
this class^. In Fig. lb the interaction potential for the 
standard XY model U (4>) = 1 — cos(0) is compared to the 
one for the Villain model at the KT-transition (T = 0.45) 
(7(0) = -Tln{^:! O oo exp(-(0-27rn) 2 /2r)}.^ 1 . The 



(a) (b) 

FIG. 2: Two groups of distinct groundstates of the 2D GF- 
FXY model, (a) When p is smaller than p c (^ 1.3479), the 
gauge- invariant phase difference <p = 7r/4 for all edges of a 
plaquette. (b) When p > p c , one edge has <j> = ir while all 
other three have = 0. The wiggled vertical lines denote the 
magnetic bond angles Ay = ir, arrows indicate phase values, 
and ± represent vortex charges. 



2D FFXY model with the Villain interaction has the 
same phase transition scenario as the usual 2D FFXY 
model i.e. a U(l) KT-transition followed by a Z2 transi- 
tion (still extremely close together but a little less close 
than for the standard 2D FFXY model) X Is this true for 
all models within the FFXY class? The answer is noi^ 3 - 
The reason is, according to us, connected to the appear- 
ance of a new groundstate. 



III. GROUNDSTATE 

Let us first consider the groundstate for the standard 
2D FFXY model on a square lattice: The spin configu- 
ration corresponding to the groundstate checkerboard is 
given in Fig. 2aJi A square with (without) a flux quanta 
is denoted by + (— )■ The arrows give the spin directions 
and the thick (thin) links are the links with (without) 
magnetic bond angles 7r (0) modulo 27r. In this configu- 
ration all the links contribute the same energy U(j) to 
the groundstate. Thus the energy for the four links con- 
stituting an elementary square is in this configuration 
4U(j). The broken symmetry of the free energy is for 
T = directly related to the fact that in order to change 
+ to — squares in Fig. 2a by continuously turning the 
spin directions from the one groundstate to the other, an 
increase of the energy is required by a finite amount for 
a number of links. This required number of links goes 
to infinity with the size of the system: the two ground- 
states are separated by an infinite energy barrier in the 
thermodynamic limit. 

The crucial point in the present context is that the 
groundstate shown in Fig. 2a does not remain the ground- 
state for all values p. As p is increased, the maximum 
link energy U(ir) decreases and at a particular value 
p c > 1 the groundstate switches to the spin configura- 
tion shown in Fig. 2b. The energy for the links around 
a square is for this configuration given by U(tt) + 3(7(0). 
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FIG. 3: Phase diagram of the 2D GFFXY model in the (p, T) 
plane. The staggered magnetization m and the helicity modu- 
lus Y give us all four combinations, all of which are realized in 
the phase diagram. The horizontal dotted line at p — 1 corre- 
sponds to the standard FFXY model which has two distinct, 
extremely close transitions. 



The critical value p c is hence given by the condition 
U(tt) + 3U(0) = 4C7(|) leading to the determination 



/ ln(3/4) 
21n(cos(7r/8) 



= 1.3479. 



(3) 



The groundstate for p > p c shown in Fig. 2b has the 
property that an infinitesimal change of the middle spin 
is enough to flip between the two checkerboard patterns 
(switching between + and — in Fig. 2b). Thus there 
is no barrier between these two checkerboard patterns 
for p > p c . This means that the broken symmetry of the 
free energy associated with the two possible checkerboard 
patterns states is restored. However, there is a new in- 
finite barrier between the two degenerate groundstates 
on opposite sides of p c : continuously turning the spins 
to change from the spin-configuration in Fig. 2a to the 
spin-configuration in Fig. 2b requires an infinite energy. 



IV. ORDER PARAMETERS 

In order to characterize the phase transition properties 
of the 2D GFFXY model one needs to identify a set of 
order parameters with which all possible transitions can 
be characterized: 

The checkerboard pattern is usually associated with 
a Z2 chirality symmetry. For T = this symmetry is 
reflected in the existence of two degenerate groundstates 
(the two checkerboards) separated by an infinite energy 
barrier. The corresponding order parameter is related to 
the staggered magnetization m defined asi^ 
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where (• • • ) is the ensemble average and the vor- 
ticity for the Ith elementary plaquette at (xi ,yi) is 
computed from s; = (1/tt) 4>ij = ±1 with 

the sum taken anti-clockwise around the given plaque- 
tte. The corresponding broken symmetry is reflected 
in the following way: for any finite system the quan- 
tity jp Y^d=\{~^-) Xl+Vl s i can with finite probability ac- 
quire any value in the range [— 1, 1] allowed by the 
model. However, in the thermodynamic limit L = 00 
only values in either the range [—1,0] or the range [0, 1] 
can be acquired. This means that the order parameter 

O = (jj^2f = i(—]-) Xl+yi si) in the tthermodynamiclimit 
only can take on the two values O = ±m. The probabil- 
ity for the two values are equal but they are separated 
by an infinite free-energy barrier. This is equivalent to 
saying that the order parameter O has a Zi symmetry 
which is broken. In the broken symmetry region m ^= 
whereas when the symmetry is unbroken m — 0. Figure 
3 shows the phase diagram in the (T,p)-plane. As seen 
m =/= corresponds to a finite region of this plane. 

For T = and p — p c , the two groundstates in Fig. 2a 
and b are degenerate and separated by an infinite energy 
barrier. For T > this should instead take the form 
of an infinite free energy barrier in the thermodynamic 
limit, separating values which a local order parameter can 
acquire. To this end one needs to identify an appropriate 
local order parameter. Such a possible order parameter 
is the defect density defined by 
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t=i 



(5) 



where the square lattice has been divided into 
—squares numerated by t where each consists of four 
elementary plaquettes. Here St is the sum of the phase 
difference around a four-plaquette s t = (l/V) Yl(ij)et $n 
which means that \s t \ can be 0,1 or 2. Thus the defect 
density can be described in the following way: Think of 
the elementary plaquettes as being either black (s = 1) 
or white (s = —1). There are always equally many black 
and white squares. The defect density measures the aver- 
age difference in the number of white and black squares 
contained in a four-plaquette. Obviously the checker- 
board groundstate corresponds to a zero defect density 
nk = 0. However, for a finite temperature the checker- 
board groundstate may contain a kink. This situation 
is illustrated in Fig. 4: start from a checkerboard pat- 
tern. The thick dotted line is a boundary between the 
two possible checkerboard patterns. The 90 degree turn 
of this line is associated with a four-plaquette with s t = 1 
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in the pressure temperature plane: the order parameter 
is the density difference on the two sides of the transi- 
tion line and the infinite free energy barrier only exists 
precisely on the transition line. 

The W(l)-symmetry is in 2D at most "quasi" broken 
because of the Mermin- Wagner Theorem™,. As a con- 
sequence the corresponding phase transitions cannot be 
described by a local order parameter. Instead the phase 
transitions can be monitored by the increase of the free 
energy caused by a uniform twist 5 of the spin angles 
across the system. Expanding the free energy F(S) for 
small values of 8 to lowest orders gives 



FIG. 4: Two checkerboard states with the boundary between 
them (denoted as thick dotted line). A kink exists where the 
boundary makes a 90-degree turn, and the kink density rik is 
measured by Eq. (5). For the four plaquettes surrounded by 
thick full line, \st \ = 1 whereas all other four plaquettes have 
even number of vortices and thus \s t \ = 0. 



which is denoted as thick solid line surrounding four pla- 
quettes in Fig. 4. Thus a kink corresponds to a defect 
with \s\ = 1 according to our definition. The defect den- 
sity defined here can be regarded as a generalization of 
the kink concept, since it does not rest on the possibil- 
ity of uniquely identifying domain boundaries. Thus the 
defect density remains a well defined concept even when 
the checkerboard symmetry is completely restored and 
to = 0. The groundstate shown in Fig. 2b is an example 
of a situation when to — 0, because switching between 
+ and - in Fig. 2b does not involve passing any energy 
barrier. Thus the defect density remains finite as T is 
lowered towards zero for any p > p c . Consequently, the 
groundstate in Fig. 2b corresponds to a finite defect den- 
sity n?j > 0. It is also obvious that the defect density is 
monotonously increasing with T. 

A phase transition associated with this order param- 
eter is signaled by either a discontinuous or a non- 
analytical behavior of rife as a function of T and p. The 
defect density makes a ddiscontinuousjump from zero to 
a finite value at p c in the limit of small temperatures and 
these two values are separated by an infinite energy bar- 
rier: the point (p, T) = (p c ,0) is the starting point of a 
phase transition line (see Fig. 3). On this phase transi- 
tion line the order parameter rik can only take on two 
values. These two values are equally probable but are 
separated by an infinite free energy barrier. Thus the or- 
der parameter rik on this phase transition line possesses 
a Zi symmetry which is broken. 

One should note that in case of n& the infinite free 
energy barrier between two different but equally proba- 
ble values of rik only resides on well defined lines in the 
(p, T)-plane, whereas the infinite barrier for the chiral- 
ity transition resides on an area of the (p, T)-plane (see 
Fig. 3). Thus the phase transition associated with the 
defect density rik is more akin to a liquid-gas transition 



S 2 S 4 
F(S)=Y-+Y 4 -. (6) 

Here, Y is the helicity modulus. It is finite in the 
low-temperature phase and zero in the high-temperature 
phased Y4 is the fourth order modulus and can be used 
to verify that the helicity modulus Y makes a discontin- 
uous jump to zero at the transition^ This discontinuous 
jump is a key characteristics of the KT-transition i 19 1 20 

V. PHASE DIAGRAM AND PHASE 
TRANSITIONS 

In Ref^ the phase transitions associated with the 
W(l)-symmetry and the ^-chirality symmetry were in- 
vestigated. The corresponding phase diagram is repro- 
duced in Fig. 3. This phase diagram has four sectors 
corresponding to all four possible combinations of transi- 
tions for a combined symmetry U(l) ® Z2: The four sec- 
tors are characterized by the four possible combinations 
(T, to) = (0, 0), (0, ^ 0),' & 0, 0), (^ 0,^0). The dashed 
horizontal line at p = 1 in Fig. 3 corresponds to the usual 
FFXY model. In this case the phase (T 7^ 0, to = 0) is 
not realized^. 

In the present paper we use all the three order pa- 
rameters described in the previous section together with 
Monte Carlo simulations in order to deduce the nature 
of the various phase boundaries. 

Fig. 5 gives a sketch of the resulting "horizontal" phase 
boundary in Fig. 3. In this blown up scale one finds that 
it has one maximum and one minimum, as well as three 
multicritical points ending three distinct phase lines. The 
critical points are denoted by A, B, and C. A fourth mul- 
ticritical point is found along the "vertical" phase line at 
much higher p and lower T (see Figs. 3 and 5). Let us first 
consider the phase boundary from T — to the critical 
point A. Across this first section of the phase boundary 
the phase transition associated with the defect density rik 
is first order. Fig. 6a illustrates the discontinuous change 
in the defect density rik- The defect-density histogram 
along this phase line has two distinct values of equal 
probability which remain distinct in the large L- limit. 
An example is given in Fig. 6b: For a given temperature 
T, the lower value corresponds to the low-p phase and the 
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FIG. 5: Magnified phase diagram near p c fa 1.3479 (com- 
pare with Fig. 5). There are in total four multicritical points 
termed A, B, C, D(see text). The critical point D shown in 
the inset occurs at much higher p and lower T 



higher to the high-p phase. As pointed out above, this 
is aanalogousto the density for a liquid-gas transition. 
Note that for T = 0.1 the p- value for the first order line 
is lower than p c (0). However, as T is further increased, 
the p-value for the first order line increases. Finally, at a 
critical temperature T c a the density difference vanishes 
with increasing system size. This is the signature of the 
critical point A which is hence the critical point ending 
the first order transition line for the defect-density. Thus 
the critical point A is analogous to the critical point end- 
ing the first order line for a gas-liquid transition. Fig. 
6c shows the defect-density histogram close to the crit- 
ical point: at the critical point the free energy barrier 
between the two phases is L-independcnt. To good ap- 
proximation this means that the ratio between the maxi- 
mum and minimum in the kink-density histogram should 
be size- independent, whereas it increases (decreases) for 
lower (higher) temperatures 2 ^. This condition is fulfilled 
to good approximation for the T-value in Fig. 6c. At 
the critical point the defect-density difference An k (the 
difference between the two maxima in Fig. 6c) should 
vanish with size as An k ~ L- f) / v 2±. The An k size scal- 
ing is shown in Fig. 6d and is consistent with an ex- 
ponent (3/v = 0.25. One can express this exponent in 
terms of the central charge c as j3/f = c/ l 22 i 23 . The 
central charge c is coupled to the symmetry of the order 
parameter. The defect density, the staggered magnetiza- 
tion and the magnetization for the 2D Ising model can all 
aacquireprecisely two distinct values with equal probabil- 
ity separated by an infinite energy barrier. The broken 
symmetry reflected by these order parameters does hence 
have a ^-character and the phase transitions are Ising 
like. The central charge is c = 1/2 for 2D Ising like tran- 
sitions. If the order parameter on the other hand is a 2D 
vector then the symmetry is £7(1) (which means that the 



order parameter with equal probability have the same 
magnitude and any direction, but that all these possibil- 
ities are separated by an infinite energy barrier) then the 
central charge is c = 1. Provided that our three order 
parameters covers all possibilities, then a phase transi- 
tion can a priori be any combination of single and joint 
transitions involving these order parameters and is hence 
contained within the designation U(l) <£> ® %2- These 
implies that the central charge can have the four values 
c = 1/2, 1, 3/2, and 2. Here a ^-transition corresponds 
to c = 1/2, an individual KT-transition or a combined 
Zi ® Z-i transition corresponds to c = 1 , the two possible 
combined U(l) <g> .^-transitions correspond to c = 3/2, 
and a combined U(\) ® Z2 ® Z2 corresponds to c = 2. 
These possibilities are tested in Fig. 6d and singles out 
c = 1 or equivalently = 1/4. This means that of 
the four possible values only c = 1 is consistent with the 
data. As will be explained below, the helicity modulus 
remains non-zero in this part of the phase diagram (com- 
pare Fig. 3) and consequently this suggests that the crit- 
ical point A reflects a combined Z2 <8> Z2 defect-density 
and chirality transition. The defect-density transition 
ends at the critical point A: as T is increased the free 
energy barrier vanishes in the large L-limit. However, 
there is a second defect-density transition line for higher 
p- values associated with aW(l) Z2 combined KT and 
defect-density transition, as illustrated in Fig. 6e and f. 
This transition is first order for higher p and ends at a 
critical point D : for a T higher than the critical point D 
there is no defect-density transition, just as for the case 
of the critical point A. 

Fig. 7a illustrates the chirality transition along the 
same phase boundary. UUp tothe critical point A (see 
Fig. 5) the transition is first order (see Fig. 7a). The 
chirality transition cannot cease at the critical point A 
because for a fixed T the free energy barrier between the 
O = ±|m| always vanishes for a large enough p. There 
are then two possibilities: it can continue alone as a Z2- 
transition or it can combine with the KT-transition into 
a joint U{\) ® ^-transition. To deduce which possibil- 
ity is the correct one, we calculate the size scaling of 
to ~ L~@/ v and decide which of the two possible sym- 
metry allowed values j3/v = c/4 = 1/8 or 3/8 is con- 
sistent with the data. Here we use standard size scal- 
ing and calculate m(T,p) for a fixed T for a sequence 
of p which crosses the phase line. As seen in Fig. 7b, a 
unique crossing point is to good approximation obtained 
for (3/v — 1/8. From this we conclude that the chirality 
transition continues alone from the critical point A as a 
.^-transition. However as we increase the temperature 
further the character of the chirality transition changes: 
using the same procedure we instead find that the value 
j3/v = 3/8 is consistent with the data (see Fig. 7c). This 
is consistent with a joint U(l) (8 Z2 KT-chirality transi- 
tion. As we increase T further we come to the critical 
point C where the KT and chirality splits up into two 
separate transitions 1 -. At this point it is possible to in- 
stead calculate the size scaling for a fixed p. The ad van- 
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FIG. 7: a) The staggered magnetization \m\ at T = 0.1. b) 
Size scaling |m| ~ L^' 3 ^" for T = 0.118 (just above point 
A), consistent with P/v = 1/8. c) Size scaling between point 
B and C, consistent with P/v = 3/8. d) Size scaling m = 
L-0/«f((T-T cC )L 1/v ) for the critical point C: Good scaling 
ccollapseobtained for T c c ~ 0.17 with v ~ 0.77 (from Fig. 
4(d)) and p/v = 3/8. 



FIG. 6: a) First order transition of n k at T — 0.1. b) Two 
valued probability distribution P(nk) at the abrupt change in 
a), c) P(rik) at the critical point A. The inset shows that the 
ratio Pmax/^min is to good approximation finite and indepen- 
dent of L. d) Finite size scaling ,Ant ~ L - ' 3 ^", is consistent 
with P/v — 1/4. e) P(rik) for a large p- value above the criti- 
cal point D. The figure illustrates that -Pmax/fmin — > oo with 
increasing system size, e) Helicity modulus transition at the 
same T and p as in e) indicating a joint first order transition. 



tage is that we can use the standard size scaling form 
m = L-P! u f{{T - T cC )L 1 /' y ). This again shows that the 
value P/v — 3/8 is consistent with the data. From this we 
deduce that there must exist a critical point B between 
A and C where the chirality transition merges with the 
KT-transition. 

Are these deductions consistent with the U{\)- 
symmetry and the helicity modulus? We argued above 
that the transition from T = to the critical point A is 
associated with the ® symmetry. This presumes 
that the W(l)-symmetry remains "quasi" broken on both 
sides of the transition, or equivalently that the helicity 
modulus T is finite on both sides. This is illustrated in 
Fig. 8a: the helicity modulus T has a minimum at the 
phase line. However, this minimum remains non-zero in 
the large L-limit, as illustrated by the inset in Fig. 8a. 
Thus T makes at most a finite jump at the transition 
and the W(l)-symmetry remains "quasi" broken. Next 
we argued that between the critical point B and C, the 
transition is a combined IA(1) ® ^-KT-chirality transi- 



tion. This means that the helicity modulus must now 
vanish at the transition. This is illustrated in Fig. 8b, 
which shows that the T-minimum now vanishes in the 
large L-limit (compare inset in Fig. 8b). Fig. 8c shows 
the same construction close to the critical point C . The 
fact that T vanishes as a power law can be verified for the 
critical point C by instead varying T for fixed p. In these 
variables the critical point C obeys a standard scaling 
relation T = L- a g(((T - T cC )lM v ) which confirms the 
power law decay of T, as opposed to the KT- universal 
jump signaling the isolated £Y(l)-transition for the XY- 
model (see Fig. 8d)^. We also note that the obtained 
critical index v « 0.77 is consistent with the data for 
to in Fig. 7d. It is also possible to use the fourth order 
helicity modulus T4 to determine the character of the 
W(l)-transition 1 ®. In Ref^ it was found from the T4- 
data that, in the interval 1.346 < p < 1.35, the character 
of the W(l)-transition was consistent with a transition 
without a discontinuous jump in the helicity modulus T. 
This is consistent with a combined U{\) ® ^-transition 
between the multicritical points B to C. 

The following picture emerges: A and D end two first 
order phase lines. A is associated with a Z2 <E> Z2- 
transition with central charge c = 1 and D with a 
U(l) <g> i^-transition with c = 3/2. B and C are both 
associated with U(l) ® ^-transitions and c = 3/2 but 
are not end-points of first order lines. 
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FIG. 8: a) The helicity modulus, T, across the phase line be- 
low point A: The inset shows that the minimum of T remains 
finite with increasing size, b) T goes to zero for large sizes 
between B and C. c) Same as b) at point C. d) Size scaling 
relation T = L" a g(((T - T cC )L 1/v ) with good ccollapse for 
v « 0.77 (a « 0.63 and T c w 0.1675 are taken from Ref-±2 



VI. STANDARD 2D FFXY MODEL 

The usual 2D FFXY model corresponds to the p = 1- 
line in Fig. 3. The critical point C for the 2D FFXY 
class is the closest multicritical point to the actual phase 
transitions of the usual 2D FFXY model (compare Fig. 
5). The critical point C is characterized by the critical 
index v » 0.77 and the central charge c = 1.5. A sin- 
gle Zi transition is characterized by v = 1 and c = 0.5. 
All the earlier papers, in which it was putatively con- 
cluded that the 2D FFXY model has a joint transition, 
the apparent value of v was in the interval 0.77 < v < 1 
(see table 1 in Refi). In particular in RefJi the val- 
ues of v and c were independently determined and given 
by v = 0.80(4) and c = 1.61(3). Thus the apparent 
multicritical point for the usual FFXY model appeared 
to have critical properties inconsistent with a single Z^- 
transition and with critical ^-values in between a single 
.^-transition and the real U(l) ® Zi multicritical point 
C for the 2D FFXY class. Furthermore, the closeness of 
the i^- values and c- values (y ~ 0.77 and c = 1.5 for C, 
respectively, v = 0.80(4) and c = 1.61(3) obtained for the 
usual FFXY model in Ref£) suggests that the putative 
multicritical point found for the 2D FFXY model is an 
artifact of the closeness to the real critical point C for 
the 2D FFXY class. 

The present consensus is that the 2D FFXY model un- 
dergoes two separate transition, a KT transition at Tkt 
followed by a ./^-transition at T% 2 with Tkt < Tz 2 X In 
particular Korshunov in Ref. 8 - has given a general argu- 
ment which purportedly states that Tkt < Tz 2 should be 



true not only for the 2D FFXY model, but also for the 
2D FFXY class studied in the present work, provided 
that the interaction is such that its groundstate is the 
broken symmetry checkerboard state. This is in contra- 
diction with the existence of the multicritical point C at 
p < p c (compare Fig. 5) which does correspond to an in- 
teraction potential with a checkerboard groundstate. We 
suggest that the reason for this fallacy of the argument is 
connected to the closeness to the (m,Y)=(0,^0)-phase. 

The most striking feature of the phase transition for 
2D FFXY model is the closeness between Tkt and Tz 2 - 
The phase diagram in Fig. 5 gives a scenario for which 
this feature becomes less surprising: The point is that 
the chirality transition and the KT-transition merge and 
cross as a function of p for the 2D GFFXY model. It 
then becomes more natural that, for some values of p, 
the transitions can be extremely close. The value p = 1, 
which corresponds to the usual FFXY model happens to 
be such a value. 

There are many other U(l) ® ^-models related to the 
2D FFXY model 1 . Although, our results only pertain to 
the 2D FFXY-class defined in this paper, we note that, 
to our knowledge, none of the phase diagrams for related 
models contain a crossing of the KT and an Ising-like 
transition. In a vast majority, the KT-transition is al- 
ways at lower temperature than the Ising-like transition 
or possibly merged. However, in the model in Refill the 
situation is reversed with the Ising-like transition below 
or merging with the KT-transition. Also in this case a 
crossing is lacking. Because there is no crossing it is no- 
toriously difficult assert whether a merging takes place or 
whether the two transitions are only extremely closed For 
example the Ising-XY model was in Refs^^ found to 
contain such a line of merged transitions. However, more 
careful MC simulations in fact suggest that the transi- 
tions are extremely close but never merge along this line. 1 
The point to note is that for our 2D GFFXY-model the 
transitions cross from which directly follows that a real 
merging exists in this case. We believe that this crossing 
is intimately related to the appearance of the additional 
groundstate. 



VII. FINAL REMARKS 

To sum up, we have found that the description of the 
phase diagram for the 2D FFXY-class of models requires 
at least three distinct order parameters consistent with 
the proposed designation U(l) <8> Z2 ® Z2: In addition 
to the usual KT U(l) transition and the chirality Z2- 
transition, there is also a defect-density transition with 
Ising like ^-character. Within our simple parametriza- 
tion of the interaction U(<j>), we have found that all 
combinations of transitions can be realized except two: 
the single ^-defect transition and the fully combined 
U(l) g> Z2 <8> ^-transition. All the others are realized, 
i.e. the single ./^-chirality transition, the single W(1)-KT 
transition, the combined ^-defect and ^-chirality tran- 
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sition, the combined Z2-chirality and U(l)-KT, and the 
combined Z 2 -defect and the U(1)-KT transition. Since 
the GFFXY-model is a subclass of the 2D FFXY-class 
this means that at least five of the symmetry allowed 
transitions can be realized. What about the remain- 
ing two? Here we speculate that a single Z2 -density 
transition will hardly be realized because it couples too 
strongly to the other transitions. However, one might 
imagine that there exists a potential U{<f) for which 
the two nearby critical points A and B are merged. 
This critical point would then correspond to a merged 
U{1) ® Z<i® .^-transition with central charge c = 2. 
We also note that Cristofano et al in Refi^, argued 



from general symmetry considerations that the full sym- 
metry of the FFXY- model allows for U{\) ® Z 2 ® Z 2 . The 
present results for the phase diagram of the 2D GFFXY 
model supports this designation. 
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